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Thermal properties of asymmetric nuclear matter are studied within a self-consistent thermal 
model using an isospin and momentum dependent interaction (MDI) constrained by the isospin 
diffusion data in heavy-ion collisions, a momentum-independent interaction (MID), and an isoscalar 
momentum-dependent interaction (eMDYI). In particular, we study the temperature dependence of 
the isospin-dependent bulk and single-particle properties, the mechanical and chemical instabilities, 
and liquid-gas phase transition in hot asymmetric nuclear matter. Our results indicate that the 
temperature dependence of the equation of state and the symmetry energy are not so sensitive to 
the momentum dependence of the interaction. The symmetry energy at fixed density is found to 
generally decrease with temperature and for the MDI interaction the decrement is essentially due 
to the potential part. It is further shown that only the low momentum part of the single-particle 
potential and the nucleon effective mass increases significantly with temperature for the momentum- 
dependent interactions. For the MDI interaction, the low momentum part of the symmetry potential 
is significantly reduced with increasing temperature. For the mechanical and chemical instabilities 
as well as the liquid-gas phase transition in hot asymmetric nuclear matter, our results indicate that 
the boundary of these instabilities and the phase-coexistence region generally shrink with increasing 
temperature and is sensitive to the density dependence of the symmetry energy and the isospin and 
momentum dependence of the nuclear interaction, especially at higher temperatures. 
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I. INTRODUCTION 

The recent progress in developing advanced radioac- 
tive beam facilities offers a great opportunity to explore 
in terrestrial laboratories properties of nuclear matter 
and/or nuclei with large isospin asymmetries. As a result, 
studies on the role of the isospin degree of freedom have 
recently attracted much attention in both nuclear physics 
and astrophysics. The ultimate goal of such studies is 
to extract information on the isospin dependence of in- 
medium nuclear effective interactions as well as the equa- 
tion of state (EOS) of isospin asymmetric nuclear matter. 
The latter, especially the nuclear symmetry energy term, 
is important for understanding not only many aspects of 
nuclear physics, but also a number of important issues 
in astrophysics [11, H, H, 13) IE B 0, Hi ■ Information about 
the symmetry energy at zero temperature is important 
for determining ground state properties of exotic nuclei 
and properties of cold neutron stars at /^-equilibrium, 
while the symmetry energy or symmetry free energy of 
hot neutron-rich matter is important for understanding 
the liquid-gas (LG) phase transition of asymmetric nu- 
clear matter, the dynamical evolution of massive stars 
and the supernova explosion mechanisms. 

Although the nuclear symmetry energy at normal nu- 
clear matter density for cold asymmetric nuclear matter 
is known to be around 30 MeV from the empirical liquid- 
drop mass formula [1, its values at other densities, 
especially at supra-normal densities, are poorly known 
P, 0. Predictions based on various many-body theo- 
ries differ widely at both low and high densities jTll, ■ 



Empirically, the incompressibility of asymmetric nuclear 
matter is essentially undetermined p^ . even though the 
incompressibility of symmetric nuclear matter at its sat- 
uration density po ~ 0.16 fm~^ has been determined to 
be 231 ± 5 MeV from nuclear giant monopole resonances 
(GMR) [13] and the EOS at densities of 2po < p < 5po 
has also been constrained by measurements of collec- 
tive flows in nucleus-nucleus collisions [3]. Fortunately, 
heavy- ion reactions, especially those induced by radioac- 
tive beams, provide a unique means to investigate the 
isospin-dependent properties of asymmetric nuclear mat- 
ter, particularly the density dependence of the nuclear 
symmetry energy. Indeed, significant progress has re- 
cently been made both experimentally and theoretically 
in extracting the information on the behaviors of nu- 
clear symmetry energy at sub-saturation density from 
the isospin diffusion data in heavy-ion collisions from 
the NSCL/MSU [H, [H, [l^. Using the isospin and 
momentum-dependent IBUU04 transport model with in- 
medium nucleon-nucleon (NN) cross sections, the isospin 
diffusion data were found to be consistent with a density- 
dependent symmetry energy of i?sym(p) ~ 31.6(p/po)''^ 
with 7 = 0.69 — 1.05 at subnormal density [ll,[i3) which 
has led to the extraction of a value of L = 88±25 MeV for 
the slope parameter of the nuclear symmetry energy at 
saturation density and a value of -ftTasy = —500 ± 50 MeV 
for the isospin-dependent part of the isobaric incompress- 
ibility of isospin asymmetric nuclear matter [l^, [l3, [3 ■ 
The extracted symmetry energy further agrees with the 
symmetry energy Egymip) = 31.6(p//Oo)°'^^ recently ob- 
tained from the isoscaling analyses of isotope ratios in 
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intermediate energy heavy ion collisions , which gives 
i w 65 MeV and i^asy ~ "453 MeV. The extracted value 
of -ft'asy — —500 ±50 MeV from the isospin diffusion data 
is also consistent with the value i^asy = — 550 ± 100 MeV 
obtained from recently measured isotopic dependence of 
the GMR in even- A Sn isotopes [l^l- These empirically 
extracted values for L and ifasy represent the best and 
most stringent phenomenological constraints available so 
far on the nuclear symmetry energy in cold asymmetric 
nuclear matte at sub- normal densities. 

Theoretically, a lot of efforts have been devoted to 
the study on the properties of cold asymmetric nuclear 
matter. However, the properties of hot asymmetric nu- 
clear matter, especially the temperature-dependence of 
the nuclear s ymm etry e nergy, h as received so far much 
less attention]2ll. [2^ . [23l [24. l25l. . For finite nuclei at 
temperatures below about 3 MeV, the shell structure and 
pairing as well as vibrations of nuclear surfaces are impor- 
tant and the symmetry energy was predicted to increase 
slightly with the increasing temperature [l^, [1^ [l^ . In- 
terestingly, an increase by only about 8% in the symme- 
try energy in the range of T from to 1 MeV was found 
to affect appreciably the physics of stellar collapse, espe- 
cially the neutralization processes . At higher temper- 
atures, one expects the symmetry energy to decrease as 
the Pauli blocking becomes less important when the nu- 
cleon Fermi surfaces become more diffused at increasingly 
higher temperatures [HI, HI [H, ^24, 25, 26]. On the other 
hand, due to the van der Waals behavior of the nucleon- 
nucleon interaction, it is expected that the so-called LG 
phase transition should occur in nuclear matter. Since 
the early work, see, e.g., Refs. [13, [Ml, [H, [s^l , many in- 
vestigations have been carried out to explore properties 
of the nuclear LG phase transition both experimentally 
and theoretically over the last three decades. For a re- 
cent review, see, e.g., Refs. 35, 36]. Most of these 
studies focused on investigating features of the LG phase 
transition in symmetric nuclear matter. New features of 
the LG phase transition in asymmetric nuclear matter 
are expected. In particular, in a two-component asym- 
metric nucleonic matter, there are two conserved charges 
of baryon number and the third component of isospin. 
The LG phase transition was suggested to be of sec- 
ond order ^37j. This suggestion together with the need to 
understand better the properties of asymmetric nuclear 
matter have stimulated a lot of work recently, sec, e.g., 
Refs. 25, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49]. 

While significant progress has been made recently, 
many interesting questions about properties of hot asym- 
metric nuclear matter remain open. Some of these ques- 
tions can be traced back to our poor understanding about 
the isovector nuclear interaction and the density depen- 
dence of the nuclear symmetry energy [l|, [3, [S] . With 
the recent progress on the constraints of the density 
dependence of nuclear symmetry energy, it is therefore 
interesting to investigate how the constrained symme- 
try energy may allow us to better understand the ther- 
mal properties of asymmetric nuclear matter. Moreover, 



both the isovector (i.e., the nuclear symmetry poten- 
tial) and isoscalar parts of the single nucleon potential 
should be momentum dependent due to the non-locality 
of nucleon-nucleon interaction and the Pauli exchange 
effects in many-fermion systems. However, effects of the 
momentum-dependent interactions on the thermal prop- 
erties of asymmetric nuclear matter have received so far 
little theoretical attention 24, 25, 26]. 

In the present work, we study systematically the ef- 
fects of isospin and momentum dependent interactions 
on the thermal properties of asymmetric nuclear matter, 
including the isospin-dependent bulk and single-particle 
properties, mechanical and chemical instability, and LG 
phase transition, within a self-consistent thermal model 
using three different interactions. The first one is the 
isospin and momentum dependent MDI interaction con- 
strained by the isospin diffusion data in heavy-ion colli- 
sions. The second one is a momentum- independent in- 
teraction (MID) which leads to a fully momentum inde- 
pendent single nucleon potential, and the third one is an 
isoscalar momentum-dependent interaction (eMDYI) in 
which the isoscalar part of the single nucleon potential is 
momentum dependent but the isovector part of the sin- 
gle nucleon potential is momentum independent by con- 
struction. We note that the MDI interaction is realistic, 
while the MID and eMDYI interactions are only used as 
references in order to explore effects of the isospin and 
momentum dependence of the nuclear interactions. 

The paper is organized as follows. In Section [Hi we 
briefly introduce the MDI, MID and eMDYI interac- 
tions and discuss some relevant thermodynamic quanti- 
ties. Results on thermal effects on the isospin-dependent 
bulk and single-particle properties of asymmetric nuclear 
matter, such as the nuclear symmetry energy, the nu- 
clear symmetry potential and isospin-splitting of nucleon 
effective mass, are presented in Section iHll The mechan- 
ical and chemical instabilities of hot neutron-rich nuclear 
matter are then discussed in Section IIVI In Section [Vl 
we present the results on the LG phase transition of hot 
asymmetric nuclear matter. All the results are calculated 
from the three interactions, and the effects of isospin and 
momentum dependence the nuclear interactions are ana- 
lyzed. Finally, a summary is given in Section [VII 



II. THEORETICAL MODELS 

A. Isospin and momentum-dependent MDI 
interaction 

The isospin- and momentum-dependent MDI interac- 
tion is based on a modified finite-range Gogny effective 
interaction [50|. In the MDI interaction, the potential 
energy density V{p, T, 6) of a thermal equilibrium asym- 
metric nuclear matter at total density p, temperature T 
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and isospin asymmetry S is expressed as follows [Ig, [50] , 
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In the mean field approximation, Eq. ([T]) leads to the 
following single particle potential for a nucleon with mo- 
mentum p and isospin t in the thermal equilibrium asym- 
metric nuclear matter, i.e., [l6l. [50| 
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Po Po 
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In the above r = 1/2 (—1/2) for neutrons (protons); 
a = 4/3; fT{r,p} is the phase space distribution func- 
tion at coordinate f and momentum p. The parameters 
Au{x),Ai{x),B,Cr,r,Cr.-T and A have been assumed 
to be temperature independent and are obtained by fit- 
ting the momentum-dependence of U{p, T — 0, S,p, r) to 
that predicted by the Gogny Hartree-Fock and/or the 
Brueckner-Hartree-Fock calculations, the zero tempera- 
ture saturation properties of symmetric nuclear matter 
and the symmetry energy of 31.6 MeV at normal nuclear 
matter density po = 0.16 fm~'^ [sOl- The incompressibil- 
ity Kq of cold symmetric nuclear matter at saturation 
density po is set to be 211 MeV. The parameters Au{x) 
and Ai (x) depend on the a; parameter according to 

2_Z? IB 

A^{x) = -95.98- a; -, Ai{x) = -120.57 -I- a; . 

cr + 1 cr -\- 1 

(3) 

The different x values in the MDI interaction are in- 
troduced to vary the density dependence of the nuclear 
symmetry energy while keepin g ot her properties of the 
nuclear equation of state fixed |l6| and they can be ad- 
justed to mimic predictions on the density dependence of 
nuclear matter symmetry energy by microscopic and/or 
phenomenological many-body theories. The last two 
terms of Eq. ([2]) contain the momentum-dependence of 
the single-particle potential. The momentum dependence 
of the symmetry potential stems from the different in- 
teraction strength parameters Ct.-t and Ct,t for a nu- 
cleon of isospin r interacting, respectively, with unlike 
and like nucleons in the background fields. More specif- 
ically, we use Cr-r = -103.4 MeV and Cr,T = -11.7 
MeV. We note that the MDI interaction has been exten- 
sively used in the transport model for studying isospin ef- 
fects in intermediate energy heavy-ion collisions induced 



by neutron-rich nuclei 3, iH, IH, [H, [H, [H [H, 113, HI . 

In particular, the isospin diffusion data from NSCL/MSU 
have constrained the value of x to be between and — 1 
for nuclear matter densities less than about 1.2po [3 [13]' 
we will thus in the present work consider the two values 
of a; = and a:: = — 1. We note that the zero-temperature 
symmetry energy for the MDI interaction with a; = and 
— 1 can be parameterized, respectively, as 31.6(p/po)*''^^ 
MeV and 31.6(p//9o)^ °^ MeV [11], and thus a; = gives 
a softer symmetry energy while x = — 1 gives a stiffer 
symmetry energy. 



B. Momentum-independent MID interaction 

In the momentum-independent MID interaction, the 
potential energy density I4iid(p, 5) of a thermally equili- 
brated asymmetric nuclear matter at total density p and 
isospin asymmetry 5 is written as 



Vmid(p, 5) = 



P P 



a p 



2 Po 1 + 7 Po'' 



+ pEryUp,^)S'- (4) 



The parameters a, (3 and 7 are determined by the in- 
compressibility Kq of cold symmetric nuclear matter at 
saturation density po [38i ] 



-29.81 - 46.90 



Ko + 44.73 
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(MeV) 
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211.05 



(5) 
(6) 
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and Kq is again set to be 211 MeV as in the MDI inter- 
action. To fit the MDI interaction at zero temperature, 
the potential part of the symmetry energy E^°l^{p, x) is 
parameterized by [IB] 

EP°Up, ^) = F{x)^ + 118.6 - F{x)] (^)«(-) (8) 
Po Po 

with F{x = 0) 129.981 MeV, G{x 0) 1.059, 
F{x = -1) ^ 3.673 MeV, and G{x = -1) = 1.569. 
We note that the MID interaction reproduce very well 
the EOS of isospin-asymmetric nuclear matter with the 
MDI interaction at zero temperature for both x — Q and 
X = —1. The single nucleon potential in the MID inter- 
action can be directly obtained as 



f/MiD(p, S, T)^a-^+ Pi-^y + U^^p, 6, r). 



Po 



Po 



(9) 



with 



u^^y{p,s,T) 



4F{x)-^ + 4(18.6 - F(a;))(-^)<=(") 



+ (18.6 - F{x))iGix) - l)(-^)^("'(52. (10) 

Po 

Therefore, the single nucleon potential in the MID inter- 
action is fully momentum-independent. It also leads to 
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the fact that the potential energy density and the single 
nucleon potential in the MID interaction arc independent 
of the temperature. 



C. Extended MDYI (eMDYI) interaction 

The momentum-dependent part in the MDI interac- 
tion is also isospin dependent while the MID interac- 
tion is fully momentum independent. In order to see 
the effect of the momentum dependence of the isovec- 
tor part of the single nucleon potential (nuclear symme- 
try potential) , we can construct an isoscalar momentum- 
dependent interaction, called extended MDYI (eMDYI) 
interaction since it has the same functional form as the 
well-known MDYI interaction for symmetric nuclear mat- 
ter JjOj. In the eMDYI interaction, the potential energy 
density Vi,MDYi(p, T, 6) of a thermally equilibrated asym- 
metric nuclear matter at total density p, temperature T 
and isospin asymmetry 5 is expressed as 



2 po 



B p 



l+a 



1 



PO J J 1 

pEP°Up,x)6'. 



p')VA2 



(11) 



Here fo{r,p) is the phase space distribution function of 
symmetric nuclear matter at total density p and temper- 
ature T. EP°l^{p, x) has the same expression as Eq. ([8]). 
We set A ^ and C = £i.^z_t£ii^ and B, a and A 

have the same values as in the MDI interaction, so that 
the eMDYI interaction gives the same EOS of asymmet- 
ric nuclear matter as the MDI interaction at zero tem- 
perature for both X = and a; = — 1. The single nucleon 
potential in the eMDYI interaction can be obtained as 



cMDYI 



(p,T,,5,p,r) = t/°(p,r,p) + C/-'^(p,<5,T), (12) 



where 



U°ip,T,p) = 



A-^ + Bi^y 

Po Po 
2C_ 

Po 



foir,P} 



(13) 



and U^^{p,5,t) is the same as Eq. (|10l) which implies 
that the symmetry potential is identical for the eMDYI 
and MID interactions. Therefore, in the eMDYI inter- 
action, the isoscalar part of the single nucleon potential 
is momentum dependent but the nuclear symmetry po- 
tential is not. For symmetric nuclear matter, the single 
nucleon potential in the eMDYI interaction is exactly the 
same as that in the MDI interaction. We note that the 
same strategy has been used to study the momentum 
dependence effects in heavy-ion collisions in a previous 
work 15211. 



D. Thermodynamic quantities of asymmetric 
nuclear matter 



At zero temperature, /t(?^, p) = W^iPfi''') ~P) ^"^^ 
the integrals in previous expressions can be calculated an- 
alytically [sot , while at a finite temperature T, the phase 
space distribution function becomes the Fermi distribu- 
tion 

frir,P)-^ (14) 



exp( — ^ J, 



1 



where fj,T is the proton or neutron chemical potential and 
can be determined from 



Pt = J fT{r,p)d^P- 



(15) 



In the above, is the proton or neutron mass and Ut 
is the proton or neutron single nucleon potential in dif- 
ferent interactions. For fixed density p, temperature T, 
and isospin asymmetry S, the chemical potential pr and 
the distribution function frir.p) can be determined nu- 
merically by a self-consistency iteration scheme [13, [fi^] . 
From the chemical potential pr and the distribution func- 
tion frif^p), the energy per nucleon E{p,T,S) can be 
obtained as 



E{p,T,6)^- 
P 



V{p, TJ) + Y,J d'p^Jrir,p} 



(16) 

Furthermore, we can obtain the entropy per nucleon 
Sr{p, T, S) as 



sAp,t,s) 



Stt 



with the occupation probability 

1 



In nT + {l—nT-) ln(l— n,-)]c?p 
(17) 

(18) 
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Finally, the pressure P{p, T, S) can be calculated from 
the thermodynamic relation 



P{p,T,S) = 



TY,Srip,T,S)-Eip,T,S) 

T 

-^PtPt- 



(19) 



III. THERMAL EFFECTS ON THE 
ISOSPIN-DEPENDENT BULK AND 
SINGLE-PARTICLE PROPERTIES OF 
ASYMMETRIC NUCLEAR MATTER 

A. Nuclear symmetry energy 

As in the case of zero temperature, phenomenologi- 
cal and microscopic studies [2l|, [S^, [23, [5^ indicate that 
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the EOS of hot asymmetric nuclear matter at density p, 
temperature T, and an isospin asymmetry S can also be 
written as a parabolic function of 6, i.e., 

E{p,T,S) = E{p,T,5 = 0) + E,yr„{p,T)S^ + O{6^). (20) 

This nice feature of the empirical parabolic law for the 
EOS of hot asymmetric nuclear matter is very useful and 
convenient for the calculation of the nuclear symmetry 
energy. With the empirical parabolic law, the tempera- 
ture and density dependent symmetry energy Egym (p, T) 
for hot asymmetric nuclear matter can thus be extracted 
from 



T. The increment of the energy per nucleon E{p,T,5) 
with the temperature reflects the thermal excitation of 
the nuclear matter due to the change of the phase-space 
distribution function /r(^, or the occupation proba- 
bility nT-(r, p). With the increment of the temperature, 
more nucleons move to higher momentum states and thus 
lead to larger internal energy per nucleon. Furthermore, 
the temperature effects are seen to be stronger at lower 
densities while they become much weaker at higher den- 
sities. At lower densities, the Fermi momentum Pfij) is 
smaller and thus temperature effects on the energy per 
nucleon E{p, T, S) are expected to be stronger. 



Esymip, T) « E{p, T,d^l)- E{p, T,S^Q). (21) 

In this sense, the symmetry energy Egymip, T) gives esti- 
mation of the energy cost to convert all protons in sym- 
metry matter to neutrons at fixed temperature T and 
density p. The parabolic approximation has been justi- 
fied for MDI interaction in Ref. [11]. We note that the 
parabolic approximation also holds very well for the MID 
and eMDYI interactions. 




FIG. 1: (Color online) Density dependence of energy per nu- 
cleon for symmetric nuclear matter (upper panels) and pure 
nuclear matter (lower panels) in MDI, MID and eMDYI in- 
teractions at r = 0, 10, 30 and 50 MeV. 

Shown in Fig. [T]is the density dependence of E{p, T, 6) 
for symmetric nuclear matter and pure neutron matter 
at T = 0, 10, 30 and 50 MeV using the MDI, MID and 
eMDYI interactions with x — and —1. For symmetric 
nuclear matter {S — 0), the parameter x — would give 
the same results as the parameter x = — 1 as we have dis- 
cussed above, and thus the curves shown in upper panels 
of Fig. [T] are the same for a; = and —1. However, 
for pure neutron matter, the parameters a; = and —1 
display different density dependence for the energy per 
nucleon E{p,T,5), which just reflects that the parame- 
ters X — and —1 give different density dependence of 
the nuclear symmetry energy as will be discussed in the 
following. From Fig. [1] one can see that the energy per 
nucleon E(p, T, S) increases with increasing temperature 
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FIG. 2: (Color online) Occupation probability distribution 
for symmetry nuclear matter and pure neutron matter at p = 
0.5po and T = 10 MeV in MDI, MID and eMDYI interactions 
with X — and x = ~1. The corresponding Fermi momentum 
Pfir) at zero temperature is indicated. 

In order to understand better the above results, we 
show in Fig. [5] the momentum dependence of the occu- 
pation probability Ut (f, p) for symmetric nuclear matter 
and pure neutron matter at p = 0.5po and T = 10 MeV 
using the MDI, MID and eMDYI interactions with x = 
and —1. For comparison, the corresponding Fermi mo- 
mentum pj(t) at zero temperature is also indicated. The 
occupation probability nrir^p) at finite temperatures is 
self-consistently determined and from which other quan- 
tities can be calculated. For symmetric nuclear matter, 
as mentioned in Section [Til the MDI interaction has the 
same occupation probability as the eMDYI interaction, 
and also the occupation probability is independent of the 
X parameter. One can see clearly from Fig. [5] that com- 
pared with the case of zero temperature, more nucleons 
move to higher momentum states at finite temperature 
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of T = 10 MeV. In addition, the results indicate that 
for symmetric nuclear matter, the occupation probabil- 
ity distribution of MID interaction is more extended than 
those of MDI and eMDYI interactions. For pure neu- 
tron matter the result of MDI interaction is very close to 
that of the MID interaction, while the eMDYI interaction 
seems to have a steeper occupation probability distribu- 
tion. Meanwhile, the value of x parameter seems to have 
little effects on the shape of the occupation probability 
distribution even though for pure neutron matter. Al- 
though here we only show the case of p = O.Spo, we note 
that these properties hold for all the densities. 

In addition, it is seen from Fig. [T]that the MDI, MID 
and eMDYI interactions give almost the same EOS for 
cold nuclear matter even at high densities though their 
parameters are constrained only at saturation density as 
discussed in Section|IIl Meanwhile, it is interesting to see 
that the MDI, MID and eMDYI interactions also produce 
quite similar EOS at finite temperatures, even at high 
densities, for both x = and x — —1. This is due to the 
fact that only low momentum parts of the single-particle 
potential is lifted at finite temperatures as we will show 
in the following while these momentum parts do not con- 
tribute much to the total energy of system through in- 
tegration with respect to the momentum. Although the 
three interactions have different occupation probability 
distributions at finite temperatures, the potential energy 
makes self-consistently the balance and leads to that they 
have the similar EOS for symmetric nuclear matter and 
pure neutron matter at finite temperatures. This feature 
implies that the interaction of the simple momentum- 
independent MID type is enough to describe the EOS of 
asymmetric nuclear matter at finite temperatures. How- 
ever, as will be discussed in following, we note that the 
chemical potential, which is a quite important quantity 
in the study of the mechanical and chemical instabilities 
as well as liquid-gas phase transition in hot asymmetric 
nuclear matter, can be significantly different for the three 
kinds of interactions. 

Now let's see the temperature dependence of the nu- 
clear symmetry energy. In Fig. [3] we show the density 
dependence of the nuclear symmetry energy at T = 0, 
10, 30 and 50 MeV using the MDI, MID and eMDYI in- 
teractions with X = and —1. For different choice of 
the parameter x = and —1, Esymip-T) display differ- 
ent density dependence with x — (—1) giving larger 
(smaller) values for the symmetry energy at lower den- 
sities while smaller (larger) ones at higher densities for 
a fixed temperature. For all the three interactions with 
both X — Q and —1, it is seen that the symmetry energy 
decreases with increasing temperature. At higher tem- 
peratures, one expects the symmetry energy Esym{p,T) 
to decrease as the Pauli blocking (a pure quantum ef- 
fect) becomes less important when the nucleon Fermi 
surfaces become more diffused at increasingly higher tem- 
peratures [n m, m, [Hi . 

Within the present self-consistent thermal model, be- 
cause the single particle potential is isospin and momen- 
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FIG. 3: (Color online) Density dependence of nuclear symme- 
try energy in MDI, MID and eMDYI interactions with x = Q 
and s = -1 at T = 0, 10, 30 and 50 MeV. 
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FIG. 4: (Color online) Temperature dependence of the total 
symmetry energy and its kinetic and potential contributions 
in MDI, MID and eMDYI interactions with x = at p = 
O.lpo, 0.5po and l.Opo. 



tum dependent with the MDI interaction, the potential 
part of the symmetry energy with the MDI interaction 
is thus temperature dependent as shown in Eq. ([1]). On 
the other hand, the potential part of the symmetry en- 
ergy with the MID and eMDYI interactions does not de- 
pend on the temperature by the construction as seen in 
Eq. (HI and Eq. (|lip . It is thus interesting to study 
how the potential and kinetic parts of the symmetry 
energy Esym{p,T) may vary respectively with temper- 
ature, which will reflect the effects of isospin and mo- 
mentum dependence of the nuclear interaction. Fig. 3] 
displays the temperature dependence of the symmetry 
energy Esym {p, T) as well as its potential and kinetic en- 
ergy parts using the MDI, MID and eMDYI interactions 
with X = at p = l.Opo, 0.5po, and O.lpo. With the 
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parameter x = —1, we note the same conclusion is ob- 
tained. For the MDI interaction, it is seen that both 
the total symmetry energy Esym {p, T) and its poten- 
tial energy part decrease with increasing temperature at 
all three densities considered. Meanwhile, one can see 
that the kinetic contribution increases slightly with in- 
creasing temperature at low temperature and then de- 
creases with increasing temperature at high temperature 
for p = l.Opo and 0.5po, while it decreases monotoni- 
cally for p = O.lpo- These features observed for the MDI 
interaction are uniquely determined by the isospin and 
momentum dependence in the MDI interaction within 
the present self-consistent thermal model. On the other 
hand, for MID and eMDYI interactions the kinetic part 
of the total symmetry energy decreases with increasing 
temperature at all the densities while the potential con- 
tribution is independent of temperature and it has the 
same value for the MID and eMDYI interactions. These 
features indicate that the temperature dependence of the 
total symmetry energy is due to both the potential con- 
tribution and kinetic contribution for MDI interaction, 
but it is only due to the kinetic contribution for the MID 
and eMDYI interactions. 

It should be mentioned that for the MDI interaction, 
the decrement of the kinetic energy part of the symmetry 
energy with temperature at very low densities is consis- 
tent with predictions of the free Fermi gas model at high 
temperatures and/or very low densities [23l . [2^ . l4ll . [42[. 
Interestingly, we can see that the temperature depen- 
dence of the total symmetry energy Egym (p, T) is quite 
similar for all the three interactions except that the MDI 
interaction exhibits a little stronger temperature depen- 
dence at higher temperatures. This is due to the fact 
that the phase space distribution function will vary self- 
consistently according to whether the single particle po- 
tential is momentum dependent or not. In addition, we 
note that as shown in Ref. [Ill, both the potential and 
kinetic parts of the symmetry energy E^ym if, T) can de- 
crease with temperature for all the densities considered 
there by using the isospin- and momentum-dependent 
BGBD interaction developed by Bombaci [ll| based on 
the well known Gale-Bertsch-Das Gupta formalism [6l| . 
The different temperature dependence of the potential 
and kinetic parts of the symmetry energy between the 
MDI and BGBD interaction is due to the fact that the 
MDI and BGBD interactions have different forms of 
the energy density functional and the MDI interaction 
leads to a more complicated momentum dependence of 
the single-particle potential. This feature implies that 
the temperature dependence of the potential and kinetic 
parts of the symmetry energy depends on the isospin and 
momentum dependence of the nuclear interactions. 



B. Nuclear symmetry potential 

To understand more clearly the effects of the isospin 
and momentum dependence of nuclear interactions, we 



discuss the single particle potential and its temperature 
dependence with the three models. In particular, we 
study the thermal effects on the nuclear symmetry poten- 
tial. The nuclear symmetry potential refers to the isovec- 
tor part of the nucleon mean-field potential in isospin 
asymmetric nuclear matter. Besides the nuclear density, 
the symmetry potential of a nucleon in nuclear matter 
also depends on the momentum or energy of the nucleon. 
In hot asymmetric nuclear matter, the symmetry poten- 
tial of a nucleon can also depend on the temperature. The 
nuclear symmetry potential is different from the nuclear 
symmetry energy as the latter involves the integration of 
the isospin-dependent mean-field potential of a nucleon 
over its momentum. Both the nuclear symmetry poten- 
tial and the nuclear symmetry energy are important for 
understanding many physics questions in nuclear physics 
and astrophysics. Various microscopic and phenomeno- 
logical models have been used to study the symmetry 
potential 7^0,11,^, 62, 63, 64, 65, 66,^, 10, 
ItiI . [zl, [tJ 74 1, and the predicted results vary widely as 



in the case of the nuclear symmetry energy. In particu- 
lar, whereas most models predict a decreasing symmetry 
potential with increasing nucleon momentum albeit at 
different rates, a few nuclear effective interactions used 
in some models give an opposite behavior. However, all 
the above studies on the nuclear symmetry potential are 
for zero-temperature and the temperature dependence of 
the nuclear symmetry potential has received so far little 
theoretical attention. The density, temperature and mo- 
mentum dependent nuclear symmetry potential can be 
evaluated from 



Usym{p,P,T) 



U„ip,p,T)-Upip,p,T) 
26 



(22) 



where Un{p,p,T) and Up{p,p,T) represent, respectively, 
the neutron and proton single-particle potentials in hot 
asymmetric nuclear matter. In the present work, we use 
S = 0.5 in the calculation for the symmetry potential 
and we note that the result can vary within only several 
percents at most by choosing different values of 6. 

In order to see the temperature effects on the nuclear 
symmetry potential, we first study the temperature de- 
pendence of the nucleon single-particle potential in hot 
nuclear matter. In Fig. we show the momentum de- 
pendence of the single particle potential of symmetric 
nuclear matter at T = 0, 10, 30 and 50 MeV and p = 
O.bpo, l.Opo and 1.5po for the MDI or eMDYI interaction. 
As mentioned in Section|TTl for symmetric nuclear matter 
the MDI and eMDYI interactions are exactly the same 
while the MID interaction is completely momentum- 
independent and thus temperature-independent as well. 
Therefore, we do not discuss here the single particle po- 
tential for the MID interaction. From Fig. [SJ it is seen 
that the single particle potentials increase with increasing 
momentum and saturate at high momenta. In addition, 
the shape of single-particle potentials are steeper with in- 
creasing density and so the momentum dependence be- 
comes stronger at higher densities. It is interesting to 



8 



see that only low momentum parts of the potentials are 
lifted with increasing temperature, and this feature is 
quite reasonable and expected since the nucleons with 
high momenta can hardly be affected by the tempera- 
ture. 
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FIG. 6: (Color online) Momentum dependence of the single 
particle potential of protons and neutrons in asymmetric nu- 
clear matter with isospin asymmetry (5 = 0.5 at p = 0.5po, 
l.Opo and 1.5po and T = 0, 10, 30 and 50 MeV in the MDI 
and eMDYI interactions with 1 = 0. 



FIG. 5: (Color online) Momentum dependence of the single 
particle potential in symmetric nuclear matter at p = 0.5po, 
l.Opo and 1.5po and T = 0, 10, 30 and 50 MeV in the MDI or 
eMDYI interaction. 

Shown in Fig. [6] is the momentum dependence of single 
particle potentials of protons and neutrons in asymmet- 
ric nuclear matter with the isospin asymmetry oi 5 = 0.5 
at T = 0, 10, 30 and 50 McV and p = 0.5po, l-Opo and 
1.5/50 for MDI (upper panels) and eMDYI (lower panels) 
interactions with a; = 0. The similar results are shown in 
Fig. [7] for X — —1. The temperature and density effects 
are seen to be very similar to that in the case of sym- 
metric nuclear matter as shown in Fig. [5l i.e., only low 
momentum parts of the potentials are lifted with increas- 
ing temperature. In contrast to the results in the case 
of symmetric nuclear matter, however, the neutron and 
proton single particle potentials in asymmetric nuclear 
matter at a fixed temperature becomes different from 
each other. For the eMDYI interaction, the potentials at 
a fixed temperature are just momentum-independently 
shifted for protons and neutrons with a higher poten- 
tial for neutrons while lower potential for protons, and 
the shifted value is sensitive to the density and the x 
value. For the MDI interaction, the isospin and momen- 
tum dependence of single particle potentials is somewhat 
complicated. At a fixed temperature, for MDI interac- 
tion with X = 0, the potential of neutrons is larger than 
that of protons at low momenta but the opposite result 
is obtained at high momenta, which indicates a steeper 
potential (stronger momentum dependence) for protons 
than for neutrons. For different x values, the single par- 
ticle potentials of MDI interaction are also seen to be 
shifted and the shifted value depending on the density 
and this is due to the fact that the term with x in Eq. ( 
[5]) is momentum-independent and depends only on the 
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FIG. 7: (Color online) Same as Fig. |6]but for x — —1. 

The above discussions imply that the nuclear symme- 
try potential for the eMDYI interaction does not depend 
on the momentum while that for the MDI interaction 
docs. Show in Fig. [5] is the momentum dependence of 
the nuclear symmetry potential at T = 0, 10, 30 and 50 
MeV and p — O.lpo, 0.5po, 1-Opo and 1.5po using the 
MDI interaction with x = and x = —\. It is seen 
that the symmetry potential decreases with increasing 
momentum for both x = and x = —\. Empirically, 
for the momentum dependence of the nuclear symmetry 
potential at zero temperature, a systematic analysis of a 
large number of nucleon-nucleus scattering experiments 
and (p,n) charge-exchange reactions at beam energies up 
to about 100 MeV has shown that the data can be very 
well described by the parametrization J/sjon = a — i'£'kin 
with a « 22 - 34 MeV and & w 0.1 - 0.2 0, [tI, S |3 . 
Although the uncertainties in both parameters a and h 
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are large, the nuclear symmetry potential at saturation 
density, i.e., the Lane potential C/Lane [Z^j clearly de- 
creases approximately linearly with increasing beam en- 
ergy. This provides a stringent constraint on the low 
energy behavior of the nuclear symmetry potential at 
saturation density. We note that for the MDI interac- 
tion with both a; = and a; = — 1 at saturation density, 
the symmetry potential agrees very well with the empiri- 
cal Lane potential. In addition, the high energy behavior 
of the nuclear symmetry potential from the MDI interac- 
tion with X — has been also found to be consistent with 
the results from recent study [13, HH in the relativistic 
impulse {t-p) approximation based on the empirical NN 
scattering amplitude [s^ or the Love-Franey NN scat- 
tering amplitude developed by Murdock and Horowitz 
[H, [8J|. On the other hand, the low energy behavior of 
the nuclear symmetry potential at densities away from 
normal nuclear density is presently not known empiri- 
cally. Experimental determination of both the density 
and momentum dependence of the nuclear symmetry po- 
tential is thus of great interest, and heavy- ion reactions 
with radioactive beams provide a unique tool to extract 
this information in terrestrial laboratories [s!], [tI, [ssj . 
From Fig. [U one can see clearly that the symmetry po- 
tentials decrease with increasing temperature, especially 
at low momenta. At high momenta (above about 500 
MeV/c), the temperature effect on the symmetry poten- 
tial is quite weak due to the fact that the nucleons with 
high momenta can hardly be affected by the temperature 
as mentioned above. 



terizes the momentum dependence of the single-particle 
potential of a nucleon. The nucleon effective mass m* is 
defined by 
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FIG. 8: (Color online) Momentum dependence of the symme- 
try potential at p = O.lpo, 0.5po, l.Opo and 1.5po and T = 0, 
10, 30 and 50 MeV in the MDI interaction with a; = and 
X = -1. 



C. Isospin-splitting of nucleon effective mass 

One of the important single-particle properties of nu- 
clear matter is the nucleon effective mass which charac- 



TJlr 



P 



rrir dUr 
p dp 



(23) 



where represents the single-particle energy while Ut 
is the single-particle potential. In such a way, the nu- 
cleon effective mass m* is related to the density of states 
m* / (2Trh)^ in asymmetric nuclear matter. By definition, 
the nucleon effective mass generally depends on the den- 
sity, isospin asymmetry of the medium, and the momen- 
tum of the nucleon [6a. Is^. Ist} . In hot nuclear medium, 
it depends on the temperature as well. At zero tempera- 
ture, when the nucleon effective mass is evaluated at the 
Fermi momentum pr = pf{T), Eq. (|23p yields the Lan- 
dau mass which is related to the /i Landau parameter of 
a Fermi liquid [H, HI, [s^l . A detailed discussion about 
different kinds of effective masses can be found in Ref. 

With the single-particle potential in Eq. ([2]), since 
the momentum-dependent part of the nuclear potential 
is independent of the parameter x, the nucleon effective 
masses are independent of the a; parameter too. The 
neutron and proton effective masses are usually different 
in asymmetric nuclear matter due to the fact that the 
momentum dependence of the single-particle potential is 
different for neutrons and protons in asymmetric nuclear 
matter. The isospin-splitting of nucleon effective mass in 
asymmetric nuclear matter, i.e., the difference between 
the neutron and proton effective masses is currently not 
known empirically [s^. We note that theoretically the 
neutron-proton effective mass splitting is still highly con- 
troversial within different approaches and /or using dif- 
ferent nuclear effective interactions [H, [tiI, [tI, [tI] . Be- 
ing phenomenological and non-relativistic in nature, the 
neutron-proton effective mass splitting in the MDI inter- 
action is consistent with predictions of all non-relativistic 
microscopic models, see, e.g., [G^ItO, Hsl; the non- 
relativistic limit of microscopic relativistic many-body 
theories, see, e.g., Refs. [H, lil, El, [El] . Recent trans- 
port model studies indicate that the neutron/proton ra- 
tio at high transverse momenta and/or rapidities is a 
potentially useful probe of the neutron-proton effective 
mass splitting in neutron-rich matter [5ll . . Since the 
momentum dependence of the single-particle potential is 
usually temperature dependent, it is thus interesting to 
see the temperature effect on the nucleon effective mass. 

In the upper panels of Fig. [9l we show the momentum 
dependence of nucleon effective mass in symmetric nu- 
clear matter at p = 0.5po, 1-Opo and 1.5po and T — 0, 
10, 30 and 50 MeV for the MDI or eMDYI interaction. 
The similar results for the neutron and proton effective 
masses in neutron-rich nuclear matter with isospin asym- 
metry (5 = 0.5 are displayed in the middle panels of Fig.[S] 
with the MDI interaction. It should be noted that for the 
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FIG. 9: (Color online) Momentum dependence of the effective 
mass of protons and neutrons in symmetric nuclear matter 
(5 = 0, upper panels) and asymmetric nuclear matter (5 = 
0.5, middle panels) at p = 0.5po, l.Opo and 1.5po and T — 
0, 10, 30 and 50 MeV in the MDI and eMDYI interactions. 
The corresponding results for the reduced isospin-splitting of 
the nucleon effective mass, i.e., {m^ — m^)/m in asymmetric 
nuclear matter [5 — 0.5) are also shown in the lower panels. 



symmetric nuclear matter, the eMDYI interaction gives 
the same nucleon effective mass as the MDI interaction. 
In addition, the effective mass is independent of the value 
of X since the the momentum dependence of the single 
particle potential does not depend on the x parameter 
as mentioned above. For the MID interaction, the single 
particle potential is completely momentum-independent 
and thus temperature- independent obviously, and the nu- 
cleon effective mass is just equal to the nucleon mass in 
free space. For a fixed temperature, one can see clearly 
from the upper and middle panels of Fig. [5] that the nu- 
cleon effective mass decreases with increasing density and 
decreasing momentum, which indicates that the momen- 
tum dependence of single-particle potential is stronger 
at higher densities and lower momenta. At high mo- 
mentum, the nucleon effective mass approaches to the 
nucleon mass in free space as the single particle potential 
is saturated. For a fixed momentum, the nucleon effec- 
tive mass increases with temperature, especially at lower 
momenta, which indicates that the temperature effect 
weakens the momentum dependence of the nuclear inter- 
action at lower momenta. In asymmetric nuclear matter 
at a fixed temperature, the neutron effective mass is seen 
to be larger than proton effective mass at a fixed mo- 
mentum, and thus leads to the isospin-splitting of the 
nucleon effective mass. 

In order to see the temperature effect on the isospin- 
splitting of the nucleon effective mass, we show the corre- 
sponding momentum dependence of the reduced isospin- 
splitting of the nucleon effective mass, i.e., (tti* —m*)/m 
in the lower panels of Fig. [SI It is indicated that the 
temperature effect on the isospin-splitting of the nucleon 
effective mass displays some complicated behaviors. At 



lower densities, the temperature effect seems to reduce 
the isospin-splitting of the nucleon effective mass for a 
fixed momentum. However, at higher densities, it de- 
pends on the momentum, i.e., at higher momenta, the 
temperature effect reduces the isospin-splitting of the nu- 
cleon effective mass while it increases the isospin-splitting 
at lower momenta. These features reflect the complexity 
of the temperature effects on the momentum dependence 
of the neutron and proton single-particle potential in hot 
asymmetric nuclear matter for the MDI interaction. It 
should be mentioned that, at zero temperature, the nu- 
cleon effective mass is usually calculated at the Fermi 
momentum Pt = pf{T), and thus the isospin-splitting is 
generally a function of the density and the isospin asym- 
metry. 



IV. MECHANICAL AND CHEMICAL 
INSTABILITIES 

The mechanical and chemical instabilities of hot asym- 
metric nuclear matter have been extensively studied 
based on various theoretical models [H [H, IH, S, [M 113, 
IHjlililill- However, effects of the momentum-dependent 
interactions on the mechanical and chemical instabilities 
have received so far little theoretical attention. In the 
following, we study the mechanical and chemical insta- 
bilities using the MDI, MID, and eMDYI interactions and 
mainly focus on the effects of the isospin and momentum 
dependence of nuclear interactions. 



A. Mechanical instability 

The mechanical stability condition for a hot asymmet- 
ric nuclear matter is 

'dP\ 

>0. (24) 

T,(5 



dp 



If the above condition is not satisfied, any growth in den- 
sity leads to the decrement of pressure. As the pressure 
is lower than its background, the nuclear matter will be 
compacted, which leads to the further growth of the den- 
sity. In such way, any small density fluctuations can grow 
and the nuclear matter becomes mechanically unstable. 

As an example to see the picture of the boundary of 
mechanical instability and the critical case, we show in 
Fig. [TO] the isothermal line in the P ^ p plane for asym- 
metric nuclear matter at T = 10 MeV and different values 
of the isospin asymmetry d using the MDI interaction 
with a; = 0. For smaller isospin asymmetries, the me- 
chanical stability condition can be violated (below the 
dotted line). When 6 is higher than the critical value 
(about 0.59 in this case), the pressure increases mono- 
tonically with density and Eq. (IM)) is satisfied for all 
densities. In the critical case, we have 



dP 
dp 



d^p 

dp^ 



= 0, 



(25) 
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FIG. 10: (Color online) Pressure as a function of density at 
fixed isospin asymmetry in the MDI interaction with x — 
at T = 10 MeV. The isothermal spinodals (ITS) and the case 
of critical isospin asymmetry are also indicated. 

where dc is the critical isospin asymmetry. Below the crit- 
ical isospin asymmetry, the extrema of the P ^ p curve 
at different isospin asymmetries form the boundary of 
the mechanical instability region, namely, the isothermal 
spinodal (ITS). 
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FIG. 11: (Color online) The boundary of mechanical (ITS) 
and chemical (DS) instabilities in the p ^ 5 plane at T = 5, 
10 and 15 MeV for the MDI, MID and eMDYI interactions 
with a; = and x — —1. 

Figs. [TT] and [T^ display the boundary of mechanical 
instabihty, i.e., ITS for the MDI, MID and eMDYI in- 
teractions at T = 5, 10 and 15 MeV with a; = and 
X ^ —1 in the p ^ 5 plane while Figs. [131 HH £^nd [15] 
display the same curves in the P ^ p plane with the in- 
clusion of the curves at constant isospin asymmetries of 
S ^ 0, 0.5, 1 and Sc- From Fig. [TT] and [T^ one sees that 
nuclear matter in the left region of the boundary of me- 
chanical instability indicated by ITS is mechanically un- 
stable, and the critical isospin asymmetry as well as the 
area of the mechanical instability region decreases with 
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FIG. 12: (Color online) Same as Fig.llllbut separately for the 
MDI, MID and eMDYI interactions to show the x dependence. 
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FIG. 13: (Color online) The boundary of mechanical (ITS) 
and chemical (DS) instabilities in the P ^ p plane at T = 5 
MeV for the MDI, MID and eMDYI interactions with x = 
and X = —1. 



increasing temperature. For each interaction, the bound- 
aries overlap at (S = for different values of x, since for 
symmetry nuclear matter the three interactions are in- 
dependent of the value of x. For the MDI and eMDYI 
interactions, the ITS has the same value at (5 = 0, as 
for symmetric nuclear matter they are exactly the same 
model as mentioned in Section JITl while for MID inter- 
action it is shifted to smaller densities at (5 = 0. The 
critical isospin asymmetry is seen to be sensitive to the 
density dependence of the symmetry energy, which can 
be seen more clearly from Figs. [13] [14] and [T5] where the 
value of 5c is indicated exactly by the dash-dotted lines. 
At T = 5 and 10 MeV the critical isospin asymmetry 
is larger for 2; = —1 than for a; = 0, while at T = 15 
MeV it is smaller for x — —1 than a; = 0. These phe- 
nomena indicate that the density dependence of nuclear 
symmetry energy and the temperature are two important 
factors to determine the critical isospin asymmetry and 
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FIG. 14: (Color online) Same as Fig.[T3]but for T = 10 MeV. 
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FIG. 15: (Color online) Same as Fig.[l3]but for T = 15 MeV. 



the area of mechanical instabihty. In addition, it is seen 
from Figs. [TT] and W2\ that the critical isospin asymmetry 
and the area of mechanical instabihty are also sensitive 
to the isospin and momentum dependence of the nuclear 
interaction, especially at higher temperatures. Detailed 
comparison indicates that the critical isospin asymmetry 
of the MDI interaction is very similar to that of the MID 
interaction at low and moderate temperatures, while it is 
similar to that of the eMDYI interaction at high temper- 
atures. Meanwhile, the area of the mechanically unstable 
region is seen to be the largest for the MID interaction 
while the smallest for the eMDYI interaction. 

In the above calculations, we fixed the temperature 
while varied the isospin asymmetry and thus got the criti- 
cal isospin asymmetry. In Fig. (THlwe fix the isospin asym- 
metry at ^ = and b — 0.5, respectively, and change the 
temperature. The case of MDI interaction with x = is 
shown as an example. It is clearly shown that increasing 
the temperature at fixed isospin asymmetry is just like 
increasing the isospin asymmetry at fixed temperature, 
and the mechanical stability condition Eq. ([M]) is satis- 
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FIG. 16: (Color online) Pressure as a function of density at 
fixed temperature in the MDI interaction with a; = at 5 = 
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temperature is also indicated. 



fied in all the densities when the temperature is larger 
than a critical temperature Tc (about 15.6 MeV at (5 = 
and 11.7 MeV at (5 = 0.5). The density at the inflection 
point, which satisfies 



dP 



0, 



(26) 



is the critical density pc, and the pressure at the inflec- 
tion point is named after as the inflection pressure Pi. 
In Fig. 1171 we show the isospin asymmetry dependence 
of the critical temperature T^, the critical density Pc, the 
inflection pressure Pi for the MDI, MID and eMDYI in- 
teractions with X = and x = —1. It is indicated that 
the critical temperature, critical density and the inflec- 
tion pressure decrease with increasing isospin asymme- 
try. Below the curves the system can be mechanically 
unstable. The critical temperature for symmetric nuclear 
matter is 15.6 MeV for the MDI and eMDYI interactions 
and 16.2 MeV for the MID interaction with both x = 
and a: = — 1. For x = Q the system is stable above a 
certain high isospin asymmetry (0.9 for MDI and MID 
model and 0.84 for eMDYI model), but for a; = -1 it 
can be mechanically unstable even for pure neutron mat- 
ter. These features indicate again that the boundary of 
mechanical instability is quit sensitive to the value of x. 
Furthermore, it is clearly seen from Fig.[T7]that the MDI 
interaction is similar to the MID interaction at low tem- 
peratures, while it is similar to the eMDYI interaction at 
high temperatures. 

To see more clearly the effect of the density dependence 
of the symmetry energy on the critical temperature, crit- 
ical density and inflection pressure, we show separately 
their isospin asymmetry dependence for each interaction 
with X = Q and a; = — 1 in Fig. [THl In each case the 
critical temperature for x — Q seems to be a little higher 
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FIG. 17: (Color online) Isospin asymmetry dependence of 
the critical temperature (upper panels), the critical density 
(middle panels), and the inflection pressure (lower panels) 
in the MDI, MID and eMDYI interactions with x = and 
X = -1. 
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FIG. 18: (Color online) Same as Fig. llTl but separately for the 
MDI, MID and eMDYI interactions with x = and x = -1. 



than that for a; = — 1 at smaller S, but vice verse at 
larger S. The critical density and the inflection pressure 
are larger for a; = than a: = — 1 at low and moderate 
isospin asymmetries, and vice verse at larger S. 



B. Chemical instability 

Now we turn to discuss the chemical instability of hot 
asymmetric nuclear matter. The chemical stability con- 
dition is 



dS 



> and 



P,T 



85 



<0, 



(27) 



and the system becomes chemical instable if either of the 
inequalities in Eq. (^7)) is violated. For example, if the 
isospin asymmetry 5 has a small growth in the region of 
chemical instability, it will even grow since more neutrons 



will move into the system from the background, because 
the low neutron chemical potential will lower the total 
energy of the whole system. This also holds true for the 
case of protons. So any isospin fluctuations will make the 
system unstable if either of the inequalities in Eq. ((27| 
is not satisfied. 




FIG. 19: (Color online) Pressure as a function of density at 
fixed isospin asymmetry in the MDI interaction with s = 
at T = 10 MeV to show how to get the chemical potential 
isobar (see text for more details). 

To analyze the chemical instability, we calculate the 
chemical potential isobar for neutrons and protons at 
fixed temperature and pressure. We do this by searching 
for the cross point between the fixed pressure line and the 
P ^ p curves with fixed isospin asymmetry. In Fig. [19] 
we can calculate the densities and the chemical poten- 
tials of the cross points, just for one isospin asymmetry at 
fixed pressure and temperature. By changing the isospin 
asymmetry from to 1 thus get the whole chemical po- 
tential isobar at a fixed temperature and pressure. We 
note that for different regions of pressure, the number 
of the cross points can be one, two or three for a fixed 
isospin asymmetry, which will be reflected in the shape 
of the resulting chemical potential isobar. The critical 
isospin asymmetry of mechanical instability is 0.59 for 
the MDI interaction at T = 10 MeV with x = and the 
corresponding curve in P ^ p plane is plotted by dotted 
line in Fig. 1191 The pressure of the inflection point is 
0.121 MeV, above which the mechanical instability dis- 
appears and only exists the chemical instability, and the 
chemical potential isobar can only be one branch for all 
values of 5. 

Fig. [20] displays the chemical potential isobar calcu- 
lated at the pressure of P = -0.05, 0.03, 0.10, 0.121, 
0.20 and 0.265 MeV/fm^. One sees that the shape of the 
curves is different for different pressures. The extrema of 
curves are just the boundary of the chemical instability, 
or diffusive spinodals (DS) , which is indicated in Fig. [^D] 
by solid circles. P = 0.265 MeV/fm'^ is the critical pres- 
sure Pn in the case of MDI interaction at P = 10 MeV 
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FIG. 20: (Color online) Chemical isobar as a function of the 
isospin asymmetry in the MDI interaction with a; = at T = 
10 MeV. The extrema of each curve are indicated. 



with a; = 0, above which the chemical potential of neu- 
trons (protons) increases (decreases) monotonically with 
b and the chemical instability disappears. The inflection 
point which satisfies 



'dS 



Pc.T 



= 



(28) 



Pc,T 



is also plotted in the figure. We note here that the ex- 
trema of fin and fip correspond to the same 6 value for 
the MDI and MID interaction and thus the critical pres- 
sure is achieved simultaneously for neutrons and protons, 
but for the eMDYI interaction the chemical potential iso- 
bar shows an asynchronous behavior for neutrons and 
protons, as will be shown in the following. This asyn- 
chronous behavior is also different for different tempera- 
tures and values of x [25^ . 

The diffusive spinodals for the MDI, MID and eMDYI 
interactions with x = and a; = — 1 at T = 5, 10 and 
15 MeV are also shown in Fig. [TT] and [12] in the p ^ 5 
plane, and in Figs. [131 [HI and [15] in the P ^ p plane. 
It is clearly shown that the diffusive spinodal envelops 
the region of mechanical instability and extends further 
out into the plane and the area of chemical instability re- 
gion decreases with increasing temperature. It should be 
noted that the region of chemical instability is between 
the curves of ITS and DS. At (5 = 0, the DS and ITS over- 
lap for the MDI and eMDYI interaction, and the curves 
with a; = and a; = — 1 overlap as well. The boundary of 
chemical instability is seen to be sensitive to the density 
dependence of the symmetry energy. At T = 5 MeV the 
maximum isospin asymmetry is larger for a; = — 1 than 
a; = 0, while at T = 10 and 15 MeV it is smaller for 
X = — 1 than a; = 0. Furthermore, it is seen from Figs.fTT] 
and [12] that the critical isospin asymmetry and the area 
of chemical instability are also sensitive to the isospin and 
momentum dependence of the nuclear interaction, espe- 
cially at higher temperatures. The maximum 5 value in 



the MDI interaction is similar to that in the MID inter- 
action at low temperature, but it becomes similar to that 
in the eMDYI interaction at high temperature. This fea- 
ture also holds in the case of mechanical instability as 
discussed before. In addition, the shape of the DS curve 
in the eMDYI interaction with a; = -latT = 5 MeV 
observed in Figs . [TTl [T^ and [T5I exhibits some unusual be- 
haviors and this is due to the asynchronous behavior of 
the chemical potential isobar for neutrons and protons. 



V. LIQUID-GAS PHASE TRANSITION 

A. Chemical potential isobar 

With the above theoretical models, we can now study 
the LG phase transition in hot asymmetric nuclear mat- 
ter. The phase coexistence is governed by the Gibbs con- 
ditions and for the asymmetric nuclear matter two-phase 
coexistence equations are 



M^(T,/,<S^) 
f^i(T,p\8^) 



P''{T,p^,S% 
M^(T,p«,J«), 
M^(T,p«,J«), 



(29) 
(30) 
(31) 



where L and G stand for liquid phase and gas phase, 
respectively. The Gibbs conditions ([29]), and for 
phase equilibrium require equal pressures and chemical 
potentials for two phases with different concentrations 
and isospin asymmetries. For a fixed pressure, the two 
solutions thus form the edges of a rectangle in the proton 
and neutron chemical potential isobars as a function of 
isospin asymmetry 6 and can be found by means of the 
geometrical construction method [H, [13, E^] ■ 

The way to calculate the chemical potential isobar is 
already explained in above, and here we take the case of 
T = 10 MeV as an example to show how we study the 
LG phase transition from the chemical potential isobar. 
The solid curves shown in Fig. [21] are the proton and 
neutron chemical potential isobars as a function of the 
isospin asymmetry 5 at a fixed temperature T = 10 MeV 
and pressure P = 0.090 McV/fm'^ by using the MDI and 
MID interactions with a; = 0. The resulting rectangles 
from the geometrical construction are also shown by dot- 
ted lines in Fig.]^!] from which one can see that different 
interactions give different shapes for the chemical poten- 
tial isobar. When the pressure increases and approaches 
the critical pressure Pq, an inflection point will appear 
for proton and neutron chemical potential isobars, i.e., 



dp 



Pc.T V / Pc,T 



= 0. 



(32) 



Above the critical pressure, the chemical potential of neu- 
trons (protons) increases (decreases) monotonically with 
S and the chemical instability disappears. In Fig. [^TJ 
we also show the chemical potential isobar at the critical 
pressure by the dashed curves. At the critical pressure. 
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FIG. 21: (Color online) The chemical potential isobar as a 
function of the isospin asymmetry 5 at T = 10 MeV in the 
MDI and MID interactions with a; = 0. The geometrical con- 
struction used to obtain the isospin asymmetries and chemical 
potentials in the two coexisting phases is also shown. 



the rectangle is degenerated to a line vertical to the S axis 
as shown by dash-dotted lines in Fig. [21] . The values of 
the critical pressure are 0.265 and 0.230 MeV/fm"^, re- 
spectively, for the MDI and MID interaction with x — 0. 
It is interesting to see that the different interactions give 
different values of the critical pressure. We note that crit- 
ical pressure is also sensitive to the density dependence 
of nuclear symmetry energy with the stiffer symmetry 
energy {x — —1) gives a smaller critical pressure [2^. 
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With the eMDYI interaction, the resulting single nu- 
cleon potential is momentum dependent but its momen- 
tum dependence is isospin- independent. Comparing the 
results with those of the MDI interaction, we can extract 
information about the effects of the momentum depen- 
dence of the symmetry potential while the effects of the 
momentum dependence of the isoscalar part of the single 
nucleon potential can be studied by comparing the results 
with those of the MID interaction. Shown in Fig. [22] is 
the chemical potential isobar as a function of the isospin 
asymmetry ^ at T = 10 MeV by using the eMDYI inter- 
action with X — 0. Compared with the results from the 
MDI and MID interactions, the main difference is that 
the left (and right) extrema of /i„ and fip do not corre- 
spond to the same 6 but they do for the MDI and MID 
interactions as shown in Fig. 1211 In particular, for the 
eMDYI interaction, the chemical potentials of protons 
and neutrons are seen to exhibit asynchronous variation 
with pressure. We note that this behavior depends on 
the temperature and the value of x. At T ~ 10 MeV 
the chemical potential of neutrons increases more rapidly 
with pressure than that of protons. At T = 5 MeV 
the chemical potential of neutrons increases more rapidly 
with a; = but vice verse with x = —1, while at T = 15 
MeV the asynchronous behavior seems not quite obvious. 
This asynchronous variation is uniquely determined by 
the special momentum dependence in the eMDYI inter- 
action within the present self-consistent thermal model. 
Actually, it is this asynchronous variation that leads to 
the fact that the left (and right) extrema of /i„ and fip 
correspond to different values of S. 

At lower pressures, for example, P = 0.090 MeV/fm^ 
as shown in Fig. [51] (a) , the rectangle can be accurately 
constructed and thus the Gibbs conditions ([55)1 , ([50)1 and 
([3T|l have two solutions. Due to the asynchronous vari- 
ation of /i„ and fip with pressure, we will get a limit- 
ing pressure above which no rectangle can be con- 
structed and the coexistence equations (|29p . ([50]) and 
([31]) have no solution. Fig. [22] (b) shows the case at 
the limiting pressure with = 0.205 MeV/fm'^ for 

X = 0. In this limit case, we note that the left edge of 
the rectangle actually corresponds to the left extremum 
of fj,p. With increasing pressure, namely, at P = 0.260, 
fin passes through an inflection point while fip still has 
a chemically unstable region and this case is shown in 
Fig. [55] (c). When the pressure is further increased to 
P = 0.305 MeV/fm^, as shown in Fig. [22] (d), fip passes 
through an inflection point while fin increases monoton- 
ically with S. We note that the asynchronous variation 
of Un and fip with pressure also depends on the value of 



B. Binodal surface 



FIG. 22: (Color online) The chemical potential isobar as a 
function of isospin asymmetry 5 at T = 10 MeV in the eMDYI 
interaction with a; = 0. 



For each interaction, the two different values of 6 cor- 
respond to two different phases with different densities 
and the lower density phase (with larger d value) defines 
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a gas phase while the higher density phase (with smaller 
S value) defines a liquid phase. Collecting all such pairs 
of 5{T,P) and S'{T,P) thus forms the binodal surface. 
Fig. [221 displays the binodal surface for the MDI, MID 
and eMDYI interactions at T = 5, 10 and 15 MeV with 
X = and a; = — 1 in the P ^ 5 plain. As expected, for 
MDI and MID interactions the binodal surface has a criti- 
cal pressure, while for the eMDYI interaction the binodal 
surface is cut off by a limit pressure. Above the criti- 
cal pressure or below the pressure of equal concentration 
(EC) point no phase-coexistence region can exist. The 
EC point indicates the special case that symmetric nu- 
clear matter with equal density coexists, which is called 
"indifferent equilibrium" [37l |. The maximal asymmetry 
(MA) also plays an important role in LG phase transi- 
tion. The left side of the binodal surface is the region of 
liquid phase and the right side the region of gas phase, 
and within the surface is the phase-coexistence region. 
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FIG. 23: (Color online) The binodal surface at T = 5, 10 and 
15 MeV in the MDI, MID and eMDYI interactions with a; = 
and X = —\. The critical pressure (CP), the limiting pressure 
(LP), the points of equal concentration (EC) and maximal 
asymmetry (MA) are also indicated. 



Now let's see what affects the critical or limit pressure 
and the region of phase-coexistence. The critical pressure 
is sensitive to the stiffness of the symmetry energy, with 
a softer symmetry energy (with x = 0) gives a higher 
critical pressure and a larger area of phase-coexistence 
region. For the case of limit pressure with eMDYI inter- 
action this holds true at T = 10 MeV and T = 15 MeV, 
but the opposite result is obtained at T = 5 MeV. For 
symmetric nuclear matter the different values of x give 
the same EC point. As the MDI and eMDYI interactions 
are the same for symmetric nuclear matter, they have the 
same EC point, but for the MID interaction the EC point 
is lower, and the amount seems to increase with increas- 
ing temperature. Below the limit pressure the binodal 
surface is quite similar for the MDI and eMDYI inter- 
actions. Comparing the MDI and MID interactions, the 
isospin and momentum dependence seems to increase the 
critical pressure in a larger amount. At T = 5 MeV and 



T = 10 MeV, the area of phase-coexistence region for 
the MDI interaction is larger than that for the MID in- 
teraction, but at T = 15 MeV the opposite result is ob- 
served. The critical or limit pressure seems not to change 
monotonically with temperature, but it is clear that the 
area of phase-coexistence region decreases with increas- 
ing temperature and the pressure of EC point increases 
with increasing temperature. 
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FIG. 24; (Color online) The phase coexistence boundary (CE) 
in P ~ p plane for the MDI, MID and eMDYI interactions 
with a; = and = — latT = 5 MeV. The isothermal 
spinodals (ITS) and diffusive spinodals (DS) are also shown 
for comparison. 
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FIG. 25: (Color online) Same as Fig.[24]but for T = 10 MeV. 

Fig. 121111111111 display the same curves (the phase coex- 
istence boundary (CE)) vn P ^ p plain, and the isother- 
mal spinodals (ITS) and diffusive spinodals (DS) are also 
included for comparison. For MDI and MID interactions 
the critical pressure is the same for chemical instability 
and binodal surface. For the eMDYI interaction as the 
binodal surface is cut off by the limit pressure, the phase- 
coexistence region can not extend the region of chemical 
instability. 
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FIG. 26: (Color online) Same as Fig. [24] but for T = 15. 



C. Maxwell construction 

The binodal surface we show in the previous section 
provide much information about the LG phase transition. 
As discussed in Ref. [s^l, we can analyze the process of 
LG phase transition in hot asymmetric nuclear matter 
by Maxwell construction. We take the case of MDI in- 
teraction with x = — latr=10 MeV as an example. 

In the left panel of Fig. [27] the system is compressed 
at a fixed total isospin asymmetry i5 = 0.5. The system 
begins from gas phase, and encountered the two-phase 
region at the point A. Then a liquid phase with higher 
density emerges from the point B, with infinitesimal pro- 
portion. As the system is compressed, the gas phase 
evolves from A to D, while the liquid phase evolves from 
B to C. In this process the gas phase and the liquid phase 
coexist and the proportion of each phase changes, but the 
total isospin asymmetry is fixed. At the point C the sys- 
tem totally changes from gas phase to liquid phase and 
leaves the phase-coexistence region. 

We can analyze this process in the phase-coexistence 
region by solving the following equations 



+ [l - X)5^ = 5p, 
A/ + (1 - A)p« - p. 



(33) 
(34) 



where S^^'^^ and p^^'^^ are the isospin asymmetry and 
density of liquid (gas) phase. The total isospin asym- 
metry 5 in this case is 0.5. The fraction of the liquid 
phase A and the total density p, from which Maxwell 
construction can be produced, can be obtained by solv- 
ing the above equations. The corresponding isotherms 
are drawn in the left panel of Fig. [2S1 The dotted line 
connecting A and C obtained by direct calculation is un- 
physical. The nearly straight line connecting A and C 
is produced by Maxwell construction and corresponds to 
the realistic process. The fraction of the liquid phase A 
from A to C is also shown in the inset and it changes 
monotonically from to 1. 
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FIG. 27: (Color online) The binodal surface at 10 MeV in the 
MDI interactions with x — —1. The points A (A') through D 
(D') denote phases participating in a phase transition. The 
critical pressure (CP), the points of equal concentration (EC) 
and maximal asymmetry (MA) are also indicated. 
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FIG. 28: (Color online) The LG phase transition process in 
the P ~ p plain for the MDI interaction with x = — 1 at 
T = 10 MeV. The system is initially compressed at fixed 
total isospin asymmetry S = 0.5 (left panel) and 5 = 0.7 
(right panel). The Maxwell construction produces the curve 
AC (A'C'). The inset displays the fraction of the liquid phase 
A from A(A') to C(C'). 



The geometry of the binodal surface offers a second 
possibility for the LG phase transition process. The sit- 
uation is displayed in the right panel of Fig. 1271 where 
the system is compressed at fixed total isospin asymme- 
try 6 — 0.7, which is larger than the isospin asymmetry 
of the CP point. The system begins from gas phase, 
and encountered the two-phase region at the point A'. 
Then a liquid phase with infinitesimal fraction emerges 
from the point B'. As the system is compressed, the 
gas phase evolves from A' to C', while the liquid phase 
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evolves from B' to D'. The system crosses the phase- 
coexistence region, but at the point C it remains at gas 
phase and leaves the binodal surface on the same branch. 
The corresponding isotherms are shown in the right panel 
of Fig.[28l The solid line rather than the dotted one con- 
necting A' and C corresponds to the real process of LG 
phase transition. In this case the fraction of the liquid 
phase A increases from to Xmax (about 0.13) and then 
drops to again. 

D. Order of LG phase transition 

In the following we consider the order of LG phase 
transition and focus on the realistic MDI interaction 
by observing the behavior of thermodynamical quanti- 
ties under fixed pressure. Here we use the pressure of 
P — 0.05 MeV/fm^ and we note that there are no quali- 
tative changes if other pressures below the critical pres- 
sure are used. We choose a relatively low pressure just 
to see more clearly the effects of the phase transition on 
the thermodynamical quantities. In the upper panels in 
Fig. [29]we show the evolution of entropy per nucleon with 
temperature at fixed pressure 0.05 MeV/fm^ for isospin 
asymmetries of <5 = and 0.5 using the MDI interaction 
with a; = and a; = — 1 respectively. The method is 
quite similar to calculating the chemical potential isobar 
above, so we do not go into the details about the calcula- 
tion method. In the upper panels in Fig. [211 the dashed 
line is obtained by direct calculation and it is unphysical, 
while the solid line corresponds to the realistic process 
and is obtained by Maxwell construction. We calculate 
the chemical potential isobar at every temperature of the 
phase-coexistence region under fixed pressure, then find 
the densities and isospin asymmetries of the coexistence 
phase from Gibbs conditions, and thus obtain the frac- 
tion of each phase by using Eqs. ([55)1 and The total 
entropy per nucleon S in the coexistence phase is then 
calculated from 

Sip,6,T) = A5^(p^<5^^) + (l-A)5'=(p'=,<S«,T), (35) 

where S''^*^'^-' can be obtained from p-^*^*^^ and (^-^^"^^ by 
using Eq. ^T7\ . From the upper panels in Fig. [21] we 
can see that a.t S — the entropy jumps at T = 10.1 
MeV, which clearly indicates that the LG phase tran- 
sition for symmetric nuclear matter under the pressure 
of 0.05MeV/fm^ (which is below the critical pressure) is 
of first order. The transition temperature is = 10.1 
MeV, and we note that this value depends on the fixed 
pressure used. On the other hand, the curves oi 5 = 0.5 
is continuous. To be more specific, we calculate the heat 
capacity per nucleon under fixed pressure from 

where the entropy per nucleon is obtained from Eq. (j35p . 
The lower panels in Fig. [^D display the heat capacity per 



nucleon as a function of temperature under fixed pressure 
0.05 MeV/fm"^ for isospin asymmetries oi S ~ and 0.5 
using the MDI interaction with a; = and a:; = — 1 respec- 
tively. For both cases of x = and x = —1 at 6 = 0.5, 
the heat capacity is continuous but its first derivative is 
not continuous, which indicates that the LG phase tran- 
sition for asymmetric nuclear matter is of second order 
according to Ehrenfest's definition of phase transitions 
[q^ . These results are consistent with those in Ref. [s^l 
with a different model. Although here we only consider 
the MDI interaction, we note that the order of LG phase 
transition does not depend on the isospin and momentum 
dependence of the nuclear interaction. 
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FIG. 29: (Color online) Evolution of entropy per nucleon and 
specific heat per nucleon with temperature under fixed pres- 
sure 0.05MeV/fm^ at 5 = and 0.5 for the MDI interaction 
with X = (left panels) and x — ~1 (right panels). In the 
upper panels, the dashed line is obtained by direct calculation 
while the solid line is from Maxwell construction. 



VI. SUMMARY 

Within a self-consistent thermal model we have studied 
in detail and systematically effects of the isospin and mo- 
mentum dependent interactions on thermal properties of 
asymmetric nuclear matter. We put the special emphasis 
on the temperature dependence of the isospin-dependent 
bulk and single-particle properties, the mechanical and 
chemical instabilities, and the LG phase transition of the 
hot asymmetric nuclear matter. In our analyses we used 
the isospin and momentum dependent MDI interaction 
constrained by the isospin diffusion data in heavy-ion 
collisions, the momentum-independent MID interaction, 
and an isoscalar momentum-dependent eMDYI interac- 
tion. Our results indicate that the EOS and the sym- 
metry energy are quite similar for the three interactions 
at finite temperature, which implies that their temper- 
ature dependence is not so sensitive to the momentum 
dependence of the interaction. In particular, the sym- 
metry energy at fixed density is found to generally de- 
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crease with temperature for the three interactions. For 
the MDI interaction, both the kinetic and potential parts 
of the symmetry energy are temperature dependent and 
the decrement of the symmetry energy with temperature 
is essentially due to the decrement of the potential energy 
part of the symmetry energy with temperature. On the 
other hand, for the MID and eMDYI interactions, the 
decrement of the symmetry energy with temperature is 
only due to the kinetic part of the symmetry energy since 
the potential part does not depend on the temperature. 
Compared with the MID and eMDYI interactions, the 
single-particle potential of the MDI interaction is isospin 
and momentum dependent, which leads to the fact that 
the symmetry potential is momentum dependent and the 
nucleon effective mass is splitted in asymmetric nuclear 
matter for the MDI interaction. It is further shown that 
only the low momentum part of the single-particle poten- 
tial and the nucleon effective mass is significantly lifted 
with increasing temperature for the MDI and eMDYI 
models. For the MDI interaction, the low momentum 
part of the symmetry potential is significantly reduced 
with increasing temperature. 

We have also analyzed the boundaries of both mechan- 
ical and chemical instabilities. We found that the area 
of both mechanical and chemical instabilities generally 
decreases with increasing temperature. Meanwhile, the 
boundaries are shown to be sensitive to the density de- 
pendence of nuclear symmetry energy. In the case of 
mechanical instability the critical isospin asymmetry is 
larger for stiffer symmetry energies (e.g., x = —1) at low 
and moderate temperatures and vice verse at high tem- 
peratures. While for the chemical instability the maxi- 
mum asymmetry is larger for stiffer symmetry energies 
at low temperatures and vice verse at high and mod- 
erate temperatures. Furthermore, it is indicated that 
the mechanical and chemical instabilities are also sensi- 
tive to the isospin and momentum dependence of the nu- 
clear interaction, especially at higher temperatures. The 
boundaries obtained with the MDI interaction are sim- 
ilar to those obtained with the MID interaction at low 
temperatures and with the eMDYI interaction at high 



temperatures. 

Finally, we have explored in detail the effects of isospin 
and momentum dependence of the nuclear interaction on 
the liquid-gas phase transition in hot asymmetric nuclear 
matter. The boundary of the phase-coexistence region 
is shown to be sensitive to the density dependence of 
the nuclear symmetry energy with a softer symmetry en- 
ergy giving a higher critical pressure and a larger area 
of phase-coexistence region. The critical pressure and 
the area of phase-coexistence region is also quite sen- 
sitive to the isospin and momentum dependence of the 
nuclear interaction by comparing the cases of MDI and 
MID interactions. For the eMDYI interaction, a limit- 
ing pressure above which the LG phase transition cannot 
take place has been found and it is shown to be sensitive 
to the stiffness of the symmetry energy as well. Fur- 
thermore, the area of phase-coexistence region decreases 
with increasing temperature, and the pressure of the EC 
point increases with increasing temperature for all of the 
three interactions considered here. The phase transition 
process and the order of phase transition have also been 
analyzed by using the Maxwell construction. Our results 
indicate that the LG phase transition for symmetric nu- 
clear matter is of first order but it becomes second order 
for the asymmetric nuclear matter. 
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